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Abstract 

The impulse approximation of pion production reactions is studied by developing a relativis- 
tic formalism, consistent with that used to define the nucleon-nucleon potential. For plane 
wave initial states we find that the usual one-body (IB) expression Oib is replaced by 02B = 
— iK(m 7T /2)OiB/rn 7T , where K(m n /2) is the sum of all irreducible contributions to nucleon-nucleon 
scattering with energy transfer of m n /2. We show that 02B ~ Oib for plane wave initial states. 
For distorted waves, we find that the usual operator is replaced with a sum of two-body opera- 
tors that are well approximated by the operator 02B- Our new formalism solves the (previously 
ignored) problem of energy transfer forbidding a one-body impulse operator. Using a purely one 
pion exchange deuteron, the net result is that the impulse amplitude for np — > d,7r° at threshold is 
enhanced by a factor of approximately two. This amplitude is added to the larger "rescattering" 
amplitude and, although experimental data remain in disagreement, the theoretical prediction of 
the threshold cross section is brought closer to (and in agreement with) the data. 

PACS numbers: 12.39.Fe, 25.40.Ve, 25.10.+S, 21.30.Fe 
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I. INTRODUCTION 



It has been known for several decades that the chiral symmetry of the strong nuclear force 
in the m q — > limit can be exploited to formulate an effective field theory using hadrons as 
fundamental degrees of freedom rather than quarks and gluons [l-3]. This theory, generically 
called chiral perturbation theory (ChPT), is widely used in both the mesonic and the A = 1 
sectors. Much effort is being put into the application of ChPT to the A = 2 sector, with 
success at low energies 4j-[6). The frontier of this program is the pion production threshold, 
where the relative momentum between colliding nucleons is p = ^Jm^m^. Pion production 
is also interesting in its own right as it provides a window into three nucleon forces ?J and 
can be used to extract information about charge symmetry breaking js, 9]. 

Being an effective theory, ChPT contains an infinite number of interactions organized in 
terms of importance according to a power counting scheme with an expansion parameter of 
m 7T /A x where A x « is the scale at which the theory ceases to become valid. For the 
problem of pion production one finds an additional parameter \ = p/rriN = wiir/p ~ 0.4. 
The fact that this parameter is large provides a significant challenge and a reorganized 
counting scheme was proposed in Ref. |l0 |. 

For a nice review of the history of meson production see Ref. 11] . The present study 
considers the specific reaction NN — > dir (the two reactions pp — > dn + and np — > dn are 
related by isospin symmetry), with the pion in an s-wave. Furthermore, we are focusing 
on the contribution of a specific diagram, the impulse approximation (IA), also known as 
"direct" production, in which the produced pion does not interact at all with the spectator 
nucleon. We would like to be clear that pion rescattering, not the IA, is known to be the 
largest contribution to the total cross section. |12J The A(1232) resonance is also known to 
contribute significantly to this observable. Our motivation for the present study is to obtain 
increased precision in the total cross section calculation and to prepare for future application 
to other observables to which the IA contributes, such as p-wave pion production. 

An additional challenge in the calculation of pion production is the presence of strongly 
interacting initial/final states. Because NN potentials are only now becoming reliable at 
such high energies, one typically employs a hybrid calculation in which a kernel is calcu- 
lated perturbatively from ChPT and then convolved with wave functions calculated from 
phenomenological potentials. Recently, this method has come under question for the IA 
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FIG. 1: Impulse approximation without initial state interactions. Solid lines represent nucleons, 
dashed lines represent pions, and double solid lines represent a deuteron. 

HQ . Ideally, one would like to derive the correct method from a relativistic formalism 
that cleanly separates effects in wave functions from those appearing in the kernel. 

Consider the IA contribution to iViV — > cLtt in the plane wave (PW) approximation where 
initial state interactions are neglected (see Fig. [1]). The amplitude for such a process has 
been estimated to go like M. lA ~ ^" cr ' Pi0G°) ~ ^\/ m N m -K ( t ) {.P)^ where <p(p) is the bound 
state wave function, evaluated in momentum space. The suppression by mj 2 was noted 
in Ref. 10[, which also included an analysis that a more detailed treatment of the power 
counting based on including initial and final state interactions introduces a power of l/m n 
via an energy denominator such that the amplitude varies as Jm^. Nevertheless, we see 
directly an explicit mj 2 times 4>(y/m m n ). 

In the physical region where m n = 140 MeV, the wave function falls as a power of 
momentum greater than unity. For small values of relative momentum, the deuteron wave 
function also falls more rapidly than an inverse power of its argument. If one takes to be 
small, the deuteron remains weakly bound 15|,ll6j and therefore its momentum wave function 
will also fall rapidly in the chiral limit. Thus the power counting can only be considered a 
very rough estimate. If we follow [loj], the impulse term is a leading order term, but the 
deuteron wave function is quite small for physical values of p and there is also a substantial 
cancellation between the deuteron s- and d-states. Thus this term's contribution to the cross 



section 12[ is small and there is a contradiction between power counting expectations and 
realistic calculations. 

This contradiction was also discussed at length in Ref. [3] where the authors intro- 
duced "wave function corrections" as a possible solution. This proposal included one-pion 
exchange (OPE) with an energy transfer of in the impulse kernel, but then subtracted 
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off a similar diagram with static OPE in order to prevent double counting. The result de- 
pended strongly on the treatment of the intermediate off-shell nucleon propagator and no 
definitive conclusion was reached. This present work is intended to settle the debate regard- 
ing the inclusion of OPE in the impulse approximation. We demonstrate, by starting from a 
consistent relativistic formalism, that non-static OPE is to be included with no subtraction 
necessary; the impulse amplitude that should be used is given in Eq. ( 1291) . Furthermore, we 
show that the traditional approach of using a one-body kernel is correct only in the absence 
of initial state interactions. 

In Sec. HIl we review aspects of the Bethe-Salpeter (BS) formalism for the two-nucleon 
problem. Section HTT1 presents the N — > Ntt operator and Sec. IIVI shows that for plane wave 
iViV — y dir, the traditional impulse approximation is approximately valid. Next, Sec. |V] 
considers the full distorted-wave amplitude by calculating the corresponding loop diagram, 
including the effects of the non-zero time components of the momenta of the exchanged 
mesons. In this section, we are able to interpret the distorted-wave amplitude as a sum 
of two-body operators. We demonstrate the new formalism by explicitly evaluating s-wave 
NN — > dir amplitudes at threshold. To aid the flow of the arguments, approximations made 
in this section are verified to be sub-leading in Appendices d| to [f|. A comparison with 
experimental cross section data is made in Sec. IVI} where we also discuss implications and 
future directions. 



II. BETHE-SALPETER BASICS 

Recall the definition of the nucleon-nucleon potential from the Bethe-Salpeter formalism. 

n 

We follow the approach of Partovi and Lomon [17[ and also consider the relationship be- 
tween the Bethe-Salpeter wave function and the usual equal time wave function as recently 



discussed in Ref. [18 ]. 



Partovi and Lomon write the Bethe-Salpeter equation for the nucleon-nucleon scattering 
amplitude A4 as 

M = K + KGM, (1) 

where K is the sum of all irreducible diagrams. The quantities Ai and K depend on the 
total four-momentum P tot and the relative four-momentum k. The two individual momenta 
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are pi 2 = -Ptot/2 ± and G is the product of two Feynman propagators: 



(■■ \ Z ~ ] \z ' -] (■■(■,. (2) 




) x - Tfijy + ie J \f 2 - m-AT + 



where is the nucleoli mass. The quantities and K differ from those of [17j by a 
factor of —i/(2n). Partovi and Lomon replace the relativistic G by the Lippmann-Schwinger 
propagator g for two particles. For scalar particles, g is obtained from G by integrating over 



the zero'th (energy) component of one of the two particles [18j|. For fermions, one must also 
project onto the positive energy sub-space of both particles. This is accomplished in the 
center of mass frame by taking 

gW P M ) = 2 n ^ - ^ • k -± ""'■■f^. - k + ^V ). (3) 

^(k)(P to t - 4m^ - 4k 2 + le) 



where E(k) = y k 2 + m^. Note that g contains the important two-nucleon unitary cut. The 
non-relativistic potential U is defined so as to reproduce the correct on-shell NN scattering 
amplitude Ai using the Lippmann-Schwinger (LS) equation 

M = U + UgM. (4) 



The quantity U is obtained by equating the M. of Eq. (CQ) with that of Eq. (j4|) to find [17] 



U = K + K{G-g)U. (5) 

In solving Eq. (TjO) for the on-energy shell scattering amplitude, U never changes the value 
of the relative energy k° away from 0. Equations (jlj) and (jSJ) are consistent with Weinberg 
power counting in which one calculates the potential using chiral perturbation theory and 
then solves the LS equation to all orders. The term G — g may be thought of a purely 
relativistic effect arising from off-shell (short-lived) intermediate nucleons, and in the present 
context a perturbative effect. 

Consider the deuteron wave function in the final state of a pion production reaction. 
For P 2 near the pole position, the second term of Eq. (TjQ) dominates and we replace the 
scattering amplitude with the vertex function V: Ai — > T, and 

T = KGY. (6) 

This equation is shown pictorially in Fig. [2j The Bethe-Salpeter wave function \1/ is defined 
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FIG. 2: Bethe-Salpeter equation near the deuteron pole. 
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FIG. 3: Bethe-Salpeter wave functions. 



(b) Deuteron 



(7) 



The wave functions of the scattering state and the deuteron are shown in Fig. |3j If one uses 
Eq. 01]), the bound-state wave function is obtained by solving the equation 



gU<j> = gUgUct>. 



(8) 



The second equation shows that U also is evaluated at vanishing values of time components 
of the relative momenta. We will treat the amplitudes ty, <f> and the state vectors \ (f>) 
(either bras or kets) as interchangeable. 

The next step is to relate \l/ with 0, which can be thought of as the usual bound-state 
wave function. This is most easily accomplished by using the projection operator P on the 
product space of two positive-energy on-mass-shell nucleons. We then have 



PG = GP = G, 



g, 



(9) 



with the last step resulting from the explicit appearance of two positive-energy projection 
operators for on-mass-shell nucleons in Eq. ([3]). We define Q = I — P and use the notation 
\&p = P\&, \1>q = and PKP = K PP , PKQ = K PQ , etc. The Q-space includes all terms 
with one or both nucleons off the mass-shell. The amplitude \l/p contains the ordinary 



nucleonic degrees of freedom so one expects that it corresponds to 0. This is now shown 
explicitly. Use / = P + Q in Eq. (JTJ) and multiply by P and then also by Q to obtain the 
coupled-channel version of the relativistic bound state equation: 

^p = GpKpp^p + GpKp Q ^ Q (10) 

V Q = GqKqpVp + GqKqqVq. (11) 

Solving Eq. ( 1TT|) for and using the result in Eq. f flUj) gives 

tf = [l-GgtfOTl-^gtfgptfp (12) 

V P = G P (Kp P + K P Q [Gq 1 - Kqq^Kqp) # p , (13) 

but one can multiply Eq. (JSJ) by P • • • P etc. to obtain the result 

U PP = K PP + KpqIGq 1 - K QQ ]- l K QP , (14) 

thus Eq. ffTBl can be re-expressed as 

■qj p = GpUpp^p = gU^p. (15) 

This last equation is identical to Eq. (JTJ). Thus we have the result that 

# P = (f). (16) 

\l/p is not the complete wave function, but we expect that is a perturbative correction 
because the deuteron is basically a non-relativistic system. 



III. THE N -> Ntt AMPLITUDE 

We now turn to the application of the Bethe-Salpeter formalism to the problem of thresh- 
old pion production. First, we remind the reader of the one-body pion production operator 



in baryon chiral perturbation theory (BChPT) 19|. For a modern review, see Ref. [20j. In 



this theory, the nucleon field is split into its heavy (H v ) and light (N v ) components, 

= e~ imNV - x {N v {x) + H v {x)) 
N v (x) = e imNV - x P + ty(x) 

H v (x) = e imNV - x P^{x) (17) 



where P± = (1 ± f)/2 and v is the velocity vector satisfying v 2 = 1 and chosen in this 
work to be v — (1, 0). The heavy component is integrated out of the path integral and the 
resulting free equation of motion for the light component has a solution, 

N(x) = ^E + m N I X I e -i(s-"w)*+i P -x ; ( lg ) 

V / 

where E = a/p 2 + m 2 N and x is a two-component Pauli spinor. In Appendix [B] we show 
that the leading order (LO) Feynman rule for the s-wave N — > Ntt amplitude vanishes at 
threshold and that the next-to-leading order (NLO) rule is 

o„ = -i^^a°C$-%T a , (19) 

Zm N Zj n 

where the derivatives act on the nucleon wave functions. 



IV. THE NN -> dir REACTION: PLANE WAVE INITIAL STATES 

Traditionally, the impulse approximation to pion production is calculated by using the 
operator of Eq. ( TT9|) as the irreducible kernel to be evaluated between non-relativistic 
nucleon-nucleon wave functions for the initial and final states. Between two-component 
nucleon spinors, 7 5 7i7° — > Oi, so 



M PW 



2m N 2f w V 



|Pi,P2>, (20) 



where the superscript on A4 indicates that we have neglected initial state interactions. Next, 
we show that Eq. (1201) is only an approximation to the full impulse amplitude derived from 
the relativistic Bethe-Salpeter formalism. We will see that this approximation is only valid 
in the absence of initial state interactions. 

For the case of plane waves in both the initial and final states, a one-body operator is 
forbidden by energy-momentum conservation, 

(P3,P4\O n \ Pl , P2 )=0, (21) 

with all the Pi on mass shell. The correct formalism must be able to explain the required 
energy transfer. Our primary thesis is that the diagram of Fig. [TJ must be obtained from 
the Feynman rules as 

= (ndabi,^) = (* \K{m n /2) G 1 0^\ Pl , P2 ) , (22) 



where G\ is the Feynman propagator of the intermediate off-shell nucleoli and K(m n /2) 
is the sum of all irreducible diagrams with energy transfer of m^/2. The second equality 
of Eq. (I22p results from the relation between T and \1/ in Eq. ([7]). This manipulation is 
necessary because (<f)\ will be used for evaluation instead of (^\, meaning that the relative 
energy must remain zero in the final state. Thus the full kernel for pion production via 
the impulse approximation is KGiO n rather than just O n . Because KGiO n is a two-body 
operator, the momentum mismatch which suppresses the IA in the traditional treatment is 
removed. 

There are two points to emphasize here. Fir stly , this treatment is not equivalent to the 



heavy meson exchange operators of Refs. 



21 



22 1 which are intended to account for the 



relativistic initial and final state interactions not present in phenomenological potentials. 
Secondly, although the assertion of Eq. f f22l greatly changes the way impulse pion produc- 
tion is calculated, one should not perform the same manipulations for the similar impulse 
approximation to photo-disintegration. The reason for this is simply that near threshold the 
nucleon remains essentially on-shell and the diagram is therefore clearly reducible. 

Next, we use \l/ = \l/p + \1/q = + \1/q and focus on the term; the other term contains 
non-nucleonic physics and may be treated as a correction. Thus the impulse approximation 
is given by 

M™ « (4>\K{mj2) G 1 0„\p 1 p 2 ). (23) 

Consider the spacetime structure of the product, GiO^- The relativistic propagator G\ 
is decomposed into three terms: 1,7°, and Y- Between two-component nucleon spinors 

7 5 7;7° -»■ &i 
7°7 5 7i7° -»■ 

f 7 5 7,7° 0, (24) 

and so we can make the replacement 

Gx O, = j ti-j + mN Qm „. £(pi) -m n + m N ^ 
(pi - q) 2 - m 2 N + ie 71 -2E(p 1 )m 7T + ml + ie 



a, (25) 



— \ 4mjy / 

where in the second line we have used that -E(pi) = tun + m n /2 at threshold. Note that 
this propagator agrees with that obtained from the Feynman rules for BChPT at LO. 
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FIG. 4: Bethe-Salpeter formalism applied to pion production for plane wave initial states. 



In order to make connection with the traditional Eq. ( TT9|) . we use the approximations 
K ~ U [corrections are 0(g — G)} and G\ ~ —i/m^ [corrections are 0(m n /mN))- Putting 
these substitutions into Eq. (|23|) . 



JV1 2B 



-a 



\viVi)- 



(26) 



The quantity U is related to the potential energy by U = — iV. Ignoring the fact that £/ 
should be evaluated for non-zero energy transfer, we use the equal-time Schrodinger equation 
to replace V — > —E d — p 2 /m N and then neglect the binding energy to find M.™ ~ -^ib^- 
This means that for a PW initial state, the traditional impulse approximation should be 
roughly adequate. This is borne out in the actual calculation of the reduced matrix elements 
for Eqs. (j50|) and (125)1. 



4 pw 

^lB 
4 PW 



-24.0 
-25.6, 



(27) 
(28) 



where we have used Ref. [Sj's definition of the reduced matrix element (we suppress the 

'or clarity) and used the same static phenomenological potential 



subscript on Ref. js|'s A 
for V (here, Argonne vl8 



23] ) that is used to calculate the wave functions. See Fig. [4] for 



a pictorial description of this section. 

It is important to note that the Bethe-Salpeter equation can also be used for the pion 



has 



rescattering diagram as shown in Fig. [5j In fact, the diagram on the right in Fig. [5] 
played an important role in the development of pion production. The authors of Ref. [24 ] 
showed that this diagram (with K approximated by OPE) becomes irreducible when the 
energy dependence of the NNirn vertex is used to cancel one of the intermediate nucleon 
propagators. This discovery resolved a problem arising from calculation of NLO loops. 
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FIG. 5: Use of the Bethe-Salpeter equation in the rescattering amplitude 
In the next section, we will show that for distorted wave initial states, Eq. ( )26|) is replaced 



by 



M DW 



2B ~ /\ 



iU(*£) iU(^) 



(29) 



where the first term contributes at leading order in the theory and the second term at 
next-to-leading order. 



V. THE NN -»■ dir REACTION: DISTORTED WAVE INITIAL STATES 
A. Definition of distorted wave operator 

There is no reason to expect the result ftA™ ~ ^° carr y over f° r a distorted 

wave (DW) initial state where p 2 = m^m^ no longer holds. Indeed, we will show that the 
traditional expression for the impulse approximation does not hold for DW amplitudes. 

The fully-relativistic initial-state wave function is denoted 



\*)i=\pi,P2) + GK\*) i 



(30) 



where the first term is exactly the initial state used in the definition of A^ PW of Eqs. ( 120]) 
and f )22|) . The complete DW impulse operator is defined as, 



M 



DW 



ISI 



(31) 



The second term includes the production operator KGiO n from Eq. ( )22|) along with initial 
state interactions, 



(32) 
(33) 
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FIG. 6: Impulse approximation using distorted waves. Solid lines represent nucleons, dashed lines 
represent pions, and ovals represent wave functions. 

where in the second line we have once again used ^ = (f) + and neglected the Q-space. 



As noted by Ref. 



251 ]. the kernel of Eq. ( 13"3"j) is a loop integral which is shown in Fig. [H] 



with K being approximated by OPE. Note that four-momenta are conserved at every vertex. 
One pion exchange is the first contribution to K in ChPT besides a short range operator 
which is irrelevant for the s-wave iViV — > dix amplitude (see Appendix |a). Nevertheless, one 
must excercise caution due to the large expansion parameter of pion production. To this 
end, we employ the deuteron of Ref. 26] which is calculated from a purely-OPE potential 
with suitable form factors. As discussed in Appendix |Cj, this deuteron wave function is 
quite accurate and increases the rescattering amplitude by only 3% over a phenomenological 
deuteron. Having then employed this deuteron wave function in the calculation of the 
traditional DW impulse approximation, we will be able to avoid any complications from 
higher order parts of the potential in our subsequent investigation of the two-body operator of 
Eq. ( J33l) . In other words, although the full potential must be present in an exact calculation, 
we expect to gain insight into the correct formalism by using an OPE-only deuteron. We 
continue to use the phenomenological potentials for the initial state. To verify that the use 
of K = OPE in the initial state does not spoil our results too much, Appendix f| examines 
heavy meson exchange in the initial state. As will be discussed, this effect is parametrically 
suppressed. 
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Note that the relative momenta of the nucleons before and after the loop (p and k) are 
external momenta to the loop integral over I = (Z°, 1), but are eventually integrated over in 
a momentum-space evaluation. Let us focus solely on the energy part of the loop integral 
and ignore the vertex factors and overall constants. We define the integral /, 

dl° 1 1 1 1 



2tt 1° - E + ie 1° 
1 



— E + ie —1° + m w — 7 — E + iel° — m n /2 — 7/2 + u>i — ie 
1 1 



(34) 



P - - 7/2 - Ui + ie 1° - + 7/2 + oj f - ie 1° - + 7/2 -u f +ie' 

where uf = (1 — p) 2 + m 2 is the on-shell energy of the initial-state pion, u'j = (1 — k) 2 +m 2 is 
the on-shell energy of the final state pion, and E = l 2 /2rriN is the kinetic energy of a single 
intermediate nucleon. Note that p 2 ~ m^rriN — •yrriN and k 2 rs —7m at. 

It is straightforward to show that if the energy components of the exchanged pions in 
the above loop are set to zero (violating conservation of four- momentum), one obtains the 
traditional impulse approximation. In this case, the pion energy denominators are pulled 
out of the integral which is then evaluated by closing the contour in the lower half plane, 



'■IB 



— uj —7 — l 2 /m N 



m^ — j — \ 2 /m N -off 



(35) 



The quantity in the first set of brackets can be recognized as the product of OPE with the 
final state wave function while the second set is the product of the initial state wave function 
with OPE. This is precisely the operator that the traditional evaluation includes. 



B. Reduction to time ordered perturbation theory (TOPT) 

Our goal is to evaluate the integral in Eq. (134)) . showing that it is a sum of TOPT terms 
which can be combined to obtain Eq. ( l29l) . To begin, we rewrite the first two factors as a 
sum, 

1 11/1 1 

1° - E + iel° - m n - E + ie ~ -m* E + ie 1° - - E + ie 

This is the key to our method because after making this split, we see two terms which each 
have the propagator structure of a rescattering box loop. Consider the first term in Eq. f l36l) ; 
this loop integral looks like a two-body operator multiplied by ^— and augmented with an 
initial-state interaction. The second term looks like the same with final-state interaction. 
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FIG. 7: Definition of the two terms in Eq. (I37p . Crosses represent the propagators which are 
absent due to the partial fractions decomposition. 
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FIG. 8: TOPT terms resulting from the I% B integral. 



We define these two integrals to be I$ B and l\ B respectively, 



'2B 



J 2B + J 2B- 



(37) 



Figure [7] illustrates the splitting described in Eq. ([3] 

Next, we perform partial fraction decomposition on each of the pion propagators, splitting 
each of the two terms into four terms. Then, we continue the decomposition process until 
each term can be expressed as a single residue. For I% B we will isolate the poles containing 
ojf and then close the contour around them (for I\ B , the Ui poles are isolated). By isolating 
the poles in this way, the resulting expression is easily recognized as the sum of six TOPT 
terms. For clarity, we show these terms pictorially for I$ B in Fig. |8] where we have left 
the overall implicit. We assume for now that the stretched box diagrams are small, as 



they were in the rescattering toy model investigation [27J and denote the sum of the four 
remaining terms with a I. 

Finally, motivated by the interpretation which is presented in the next section, we alge- 
braically re-combine these four terms to find 



ja 
J 2B 



l 2B 



1 



1 + 



1 



-oj) -IE - 7 



u f + 5 b (l) 



1 



Ui + 5 a (\) J m n - 2E - 7 -co* 



(38) 



(m 7r /2)2-(a; i + 5 6 (l))^ 



(39) 
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where we have separated out terms involving S a and S\ 



b- 



because (as will be shown in the next section) they are sub-leading and we will neglect them 
in the main body of this work. The only approximation made in the evaluation of the loop 
integral to obtain Eqs. ([38]) and ([39]) is to neglect the stretched boxes. Let us pause to 
summarize what we have done so far: (1) the DW amplitude was written down as a loop 
integral, (2) partial fractions was used to split the product of the two nucleon propagators 
into a sum J| B + I% B , (3) the loop integrals were evaluated and the result expressed in terms 
of TOPT diagrams, and (4) the TOPT diagrams were algebraically combined into a form 
useful for the following interpretation. 



C. Interpretation 

Although not obvious at first sight, convolution of the operator corresponding to Eq. ( [38]) 
with wave functions as defined in Eqs. fl3TT) and (1331) results in an amplitude approximately 



equivalent to that which one obtains by using the operator shown in Fig. 9(a) The same 



is true of Eq. fl39|) with Fig. 9(b) , and together they replace the traditional (one-body) 
impulse approximation with Eq. ( l29l) . Furthermore, the operator that results from Eq. f )39l 
is expected to be small by power counting arguments. The task of this subsection is to verify 
these statements in detail. 

In Eq. ([3*8]) the factor (m^ — 2E — is interpreted as the product of the 

two-nucleon initial-state wave function with static OPE. This is the statement that 

1 1 



gU°™. (42) 



m n — 2E — 7 — ujf 

This factor can be absorbed (after adding in the PW term) using the zero-relative-energy 
Lippmann-Schwinger equation that is employed by the phenomenological potentials we are 
using. We will continue to refer to the initial wave function as a function of p and p«, so 
absorbing this factor means that we set 1 = p. 

Likewise, in Eq. ([3*9]) the factor (— Wj) _1 (— 2E — 7)- 1 is interpreted as the product of 
static OPE with the two-nucleon final-state wave function: U OPE g. Absorbing this factor 
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FIG. 9: Two-body impulse production operators (pion exchange is non-static). 

into the wave function, we set 1 = k. The remaining factors of / 2B and / 2B become the 
two-body impulse production operators, 



U 2B 



O2B 



cti (p - k)cr 2 ■ (k - p) 

(m^/2) 2 -(u; f + cUp)) 2 

4(k) 



1 + 



Sk 



1 



u f + 5 b (k) 



1 



1 + 



Safe) 

Uf 

S b (k) 



-ra„ 



$a ( P ) 



+ S a (p) 



UJ, 



cr 1 ■ (p - k)cr 2 ■ (k 



Sp 

P) 



{mj2f - (co t + 6 b (k)y 



(43) 



(44) 



where we have now made explicit the momentum dependences of the vertices and used 
S = (<Ti + er 2 )/2. It is also important to include form factors in the OPE which match those 
of the wave functions. These form factors are present in our calculation even though we 
leave them out of this expression for the sake of generality. 

Next, note that in the evaluation of the matrix element using Eq. (|43p . the initial state 
wave function is peaked about its plane wave value p ~ p^, and thus E m n /2 — 7/2 and 
5 a (p) w 0. On the other hand, in Eq. ( jUJ) , we have k w kj and E m —7/2 and 5&(k) w 0. 
If we were to neglect all the <5's, we would have 

cr 1 ■ (p - k)cr 2 ■ (k - p) 1 



2L>> 



m n /2) 2 - ((p - k) 2 + m 2 ) -m^ 
1 o"! ■ (p - k)cr 2 • (k - p 



Sp 



Sk 



(45) 



(46) 



(m n /2) 2 - ((p - k) 2 + m 2 ) ' 

which suggests that these operators can be approximately interpreted as the diagrams in 
Fig. Thus we have finally obtained our central result [Eq. (|2"§|) ] which states that the 
correct impulse approximation is a two-body operator. The contribution to pion production 
given in Eq. ( )29|) is not replacing the rescattering diagram (which is also two-body), but 
rather replacing the traditional contribution which has been referred to as the impulse 
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TABLE I: Threshold reduced matrix elements calculated with an OPE deuteron and various phe- 
nomenological initial states. The first row shows the traditional impulse approximation (one-body) 
while the second and third show our replacement (two-body). 



Avl8 Reid '93 Nijm II 



aDW 
A 1B 


8.3 


7.1 


5.4 


iDW,a 


17.4 


13.5 


7.8 


iDW,b 
A 2B 


-1.5 


-2.2 


-6.9 



approximation (or direct production). Note that if we assign standard pion production 
power counting to these diagrams, Fig. 9(a) is O (^yj^^j while Fig. 9(b) is O (j^^j- 111 
the next section the approximate expressions given in Eqs. (|45p and fj46|) are numerically 
evaluated. Nevertheless, we acknowledge the importance of verifying that the 5 terms are 
indeed small and relegate that discussion to Appendices |D| and|E. 



D. Evaluation of two-body operators 



Next, we calculate the threshold s-wave np — > dn amplitudes corresponding to Eqs. (1431) 
and ( 1461) . We do not present the details here as most are given in Ref. 13j. Again, we remind 
the reader that for the sake of consistency we use a deuteron wave function calculated from 
a purely-OPE potential (with form factors as described in Appendix [C]). For the initial-state 
distorted waves, we use three different phenomenological potentials (Avl8 {2^, Nijmegen II 
2^ |, and Reid '93 28]). In Table HJ we display the results in terms of the reduced matrix 
elements of Ref. 

The first row of Table [J give s the traditional (one-body) impulse approximation, which 
is slightly bigger than Ref. 13j due to the use of the OPE deuteron. The next row shows 
that the new two-body operator (at leading order) is roughly twice as large as the tradi- 
tional calculation it is replacing. We mention here that the significant cancellation between 
deuteron s- and <i-states remains, keeping the impulse amplitude smaller than rescattering; 
however, the cancellation is less complete when using our new two-body operator. The final 
row verifies that the 0^ B diagram is smaller than the O^-q diagram, as dictated by the power 
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counting. The Nijmegen II potential provides a bit of deviation from these results, and it 
will be interesting to investigate other potentials to determine the true model dependence 
of this calculation. In finding these results, it is important that the pion propagators of 
Eqs. (|45p and (|46p be implemented in a manner consistent with the potential used for the 
wave function of Fig. Namely, the cutoff procedure of the convolution integral with form 
factors needs to match that by which the potential was constructed. Appendix [C] contains 
the details of this procedure. 

Our conclusion is that the traditional impulse approximation is an underestimate. While 
it is true that several approximations were made in order to permit final expressions as 
simple as Eqs. (j4"5|) and (l4"5|) . we believe this conclusion to be sound. The 5 terms do not 
defy their classification as sub-leading (see Appendices D and E), and Appendix shows 
that using K = OPE in the initial state is at least reasonable. In summary, we simply claim 
that Eq. (145]) is the new impulse approximation at leading order in the effective field theory. 
The corrections in the aforementioned appendices, in addition to Eq. ( )46|) contribute to 
the next-to-leading order calculation, which needs to be systematically considered in a later 
work. 

Finally, it is important to note that although the OPE deuteron reproduces the phe- 
nomenological results for the rescattering diagram quite well, the numbers in this section 
are greatly changed if a phenomenological deuteron is used. Using Avl8 we find A^q = 4.9, 
and by using the cutoff procedure of Avl8 for the two-body operators, we find A™' a = 33.5, 
^2B V t> = — 2.8. Thus, the ratio of our new two-body operator to the traditional impulse 
operator is ~ 7 instead of the ~ 2 presented above. At this time one is faced with a choice 
of either: (1) using a "correct" phenomenological deuteron and leaving out parts of the 
potential when calculating the two-body kernel or (2) using an inexact OPE deuteron with 
a completely self-consistent kernel. For the time being, we believe the latter to be more 
trustworthy, if not ideal. 

VI. DISCUSSION 

Experimental data for pion production near threshold are reported in terms of two pa- 
rameters, a and (3, defined for np — > dn°, 

a(ri) = \{oai + / V) , (47) 
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TABLE II: Threshold reduced matrix elements extracted from experiment 



Experiment ^expt 

np^d-K [29] 80.1 ± 1.1 

pp -> dir + (Coulomb corrected) [30] 85.2 ± 1.0 

pp — > dn + (Coulomb corrected) [31] 84.6 ± 1.9 

Pionic deuterium decay [32] 93.8+2 o 



TABLE III: Rescattering (RS) and total reduced matrix elements for a variety of potentials. The 
second line shows the traditional calculation (with a one-body IA) while the third shows our 
replacement (with a two-body IA). 

Avl8 Reid '93 Nijm II 
RS 69.8 72.1 74.0 

RS + IA (IB) 78.1 79.2 79.4 
RS + IA (2B) 87.2 85.6 81.8 



where rj is the pion momentum in units of its mass. Table VI of Ref. 13[ shows the results 
obtained by the four most recent experiments. Since the present calculation is performed at 
threshold (r) — 0), we compute only the value of a, 



a 



^zl^l 2 , (48) 



128vr 2 sp 

where s = (rrid + m^) 2 is the square of the invariant energy. For ease of comparison, we 
invert Eq. (l4~8j) . plug in the results of the mentioned experiments, and propagate the errors 
to find Table M 

The full theoretical amplitude includes not only the impulse diagram but also the rescat- 
tering diagram, which is given in Table IIHI along with the total amplitude using either the 
traditional one-body or the leading-order two-body impulse diagram. The uncertainty in an 
effective field theory calculation is estimated by the power counting scheme. In this work, we 



have included both the rescattering and the impulse diagrams up to O y^/m^/rriNj ■ There- 
fore one might assign an uncertainty of m n /m N = 14% to the calculation but stress that 
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such an estimate based solely on power counting is rough at best. Taking this uncertainty, 
we see that the theory update presented here changes the situation from under-prediction 
of the most recent pionic deuterium experiment by ~ 1.3 a, to under-prediction by ~ 0.7 a. 

In summary, we have developed a consistent formalism that allows one to separate effects 
of the kernel from those of the wave functions, finding a new impulse approximation kernel. 
This two-body operator, given in Eq. ( 129|) . replaces the traditional one-body impulse ap- 
proximation and is the central result of the present work. We numerically investigated the 
simplest example (s-wave NN — > dn) and found the impulse amplitude to be increased by a 
factor of roughly two over the tradational amplitude. This calculation was performed with a 
regulated OPE deuteron which has advantages and disadvantages as described in the body 
of this work. Rescattering remains the dominant contribution to the cross section. We find 
that the updated total cross section is ~ 10% larger than before and is in agreement with 
experiment at leading order. We verified that corrections to the new impulse approximation 
(which together with other loops and counterterms will contribute at next-to-leading order) 
do not destroy these results. 

These findings suggest several directions for future research. Firstly, one needs to develop 
a power counting scheme for the "Q space" discussed in Sec. [Ill Secondly, the significant 
model dependence of the new formulation of the impulse approximation needs to be inves- 
tigated in a renormalization group invariant way. Thirdly, it will be very interesting to see 
the impact of this increased impulse amplitude on the pp — > ppn cross section which is sup- 
pressed due to the absence of rescattering. Finally, one could look at the energy dependence 
(p-wave pions) of NN — > NNtt, for which there is an abundance of experimental data. 
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TABLE IV: Parameters used. 



m n = 134.98 MeV g A = 1-32 MeV 
m N = 938.92 MeV f n = 92 A MeV 



Appendix A: Lagrange densities 

We define the index of a Lagrange density to be 

v = d+ f --2, (Al) 

where d is the sum of the number of derivatives and powers of m w , and / is the number of 
fermion fields. This represents the standard power counting for nuclear physics. The v = 



Lagrangian (with spatial vectors in bold font) is 10] 



£ (0) = \ (d7r a ) 2 - \m 2 y a + Nkd N + |f M (r a cr ■ Vvr a ) N + (A2) 

Z Z ^Jtt 

where r a and a are the Pauli matrices acting on the isospin and spin of a single nucleon. 
The "+..." indicates that only the terms used in this calculation are shown. 

The v = 1 Lagrangian includes recoil corrections and other terms invariant under 
SU(2) L xSU(2) i? . 



|f irfraw ■ VN + H.c. 



+ ..., (A3) 



2m A? 2m n 

where we use the values given in Table IIVL Note that the terms with the q low energy 
constants which appear at this order do not get promoted in MCS for these kinematics and 
are thus not used. Also, the terms with the di low energy constants do not contribute to 
s-wave pion production. Finally, the NNNN contact terms, C$,Ti do not contribute because 
we are using a potential with a repulsive core [Ri{r)Rf{r) — > as r — > for li — 1, If — 0]. 

Appendix B: iV — )• Ntt from BChPT 

The LO NNir interaction reads, 

£(o) c n^ ± ^ Nj (Bi) 
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where u± tfM = — v ■ u v^, = i{v)d^u — ud^), and u 2 = e iTana ^ 7r . We find 



u. 



J it ^Jtt Jit 



Tn 



i = 1,2,3 



i ± = 7o(«0 - U ■ 1) - 7i (^--^^^a 
= -7-liOi-Ka, 

JtT 

and thus the Feynman rule for an outgoing pion with momentum q and isospin a is 

0^ = -i[^-lrf){iqi)T a . 



fB2l 



(B3) 



At threshold, the pion four-momentum is q = (m n , 0, 0, 0) making O[ 0) = 0. This reflects 
the fact (well-known from current algebra) that threshold pion production proceeds via the 
off-diagonal, and therefore l/m N suppressed, interaction 5' 7r 7 5 7°g°r. In the effective theory, 
this recoil correction shows up in the NLO Lagrangian 



2m N 

where the spin vector is = — -ol^il^f — v^). Thus the Feynman rule is 



(B4) 



0« = -i 



- — (-477V 
Jtt 

■ m n 9A 5 0/^ &\ 

- % 7> TT7 7i7 (V - V 

2m N 2/^ 



(B5) 



where the derivatives act on the nucleon wave functions. 



Appendix C: One pion exchange deuteron 

In this Appendix, we present the method by which the deuteron wave function is calcu- 
lated for use in Sec. [V] This method is taken directly from the work of Friar, Gibson, and 



Payne 26j. The OPE potential is defined to have central (Y) and tensor (T) parts, 



K(r) = / 2 m T ^ [a, ■ * 2 Y(r) + S 12 T(r)} , (CI) 

where f 2 = 0.079 (to be distinguished from f n ) measures the strength of the pion-nucleon 
coupling and Su is the standard tensor operator. The deuteron has isospin zero and spin 
one, so we have 

K(r) = -/ 2 77v [Y(r) + S 12 T(r)\ . (C2) 
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The Y and T functions are expressed as derivatives of the Fourier transform of the pion 
propagator, 

Y(r) = h'^x)-h' (x)/x 
T(r) = h^(x)+2h' Q (x)/x 

= whS^^h/'^ (C3 > 

where x = m^r and F^nn is the form factor for which we use, 

'A 2 -m 2 ^ n 



In Ref. 



26j . it is shown that 



Y(r) = —- P 3 e-? x V \ (5i(Px) - 2i6 i - 1 (J3x)) (C5) 

i=0 



T(r) = — f 1 + - + 1) - P 3 e-P* V £ [Wx) - (2z - ^(/Jx)] (C6) 

v 7 i=0 

where (3 = A/m^ and £ = (/3 2 — l)/2/3 2 and the Si are defined by 

S i+1 ((3x) = (2i - l)<Ji (fix) + (px)Hi- X {px) (C7) 



along with 5 = 1 / /5x and #i = 1. One of the results of Ref. 26] is that larger values of n 
lead to better fits to experimental data. We use n = 5 and /3/vTo = 5.687805 in order to 
precisely reproduce the binding energy Eg = 2.2246 MeV. The wave functions are calculated 
by integrating in from r max = 100 fm and adding together two linearly independent solutions 
such that the sum vanishes at r m i n = 0.01 fm. As shown in Fig. [TOj the results are close to 
the "correct" Avl8 deuteron. 

In Table |VJ we display the quadrupole moment and mean square charge radius of Avl8, 



this OPE potential, and experiment (as quoted in [261]) . It is clear that the form factors 
in the OPE potential make it difficult to distinguish this construction as less accurate than 
Avl8. Finally, in Table IVIj we display the reduced matrix elements for the rescattering pion 
production diagram evaluated with both the phenomenological potentials and the deuteron 
of this section. Since this diagram makes the largest contribution to the cross section we 
need to verify that neglecting non-OPE parts of the potential does not dramatically change 
this amplitude. Indeed, we observe what should be expected: since the rescattering diagram 
is not as sensitive to the core of the deuteron, using the OPE wave function in place of the 
standard one has only a small effect. 
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r(fm) 



FIG. 10: (Color online) Deuteron s- and d-state wave functions (the s-state is larger). The poten- 
tials used to calculate the wave functions are Avl8 (black) and the cutoff OPE described in this 
section (red). 



TABLE V: Deuteron properties. 



Potential Q (fm 2 ) (r 2 ) 1 / 2 (fm) 

Avl8 0.270 1.968 

OPE (n = 5) 0.282 1.939 

Experiment 0.2859(3) 1.955(5) 



Appendix D: Effect of the 5 terms: diagram 

In this section we calculate the correction terms to the first two-body DW amplitude [Eq. 



13]) ] which is shown in Fig. 9(a) Assuming that the 5's truly are small compared to Eq. 
(145]) . we will only worry about calculating them one at a time, numbering the contribution 
of the S's from right to left as 1, 2 and 3. Note that we will display the results as calculated 



TABLE VI: Effect of using OPE deuteron on rescattering diagram. 



Deuteron 


Avl8 


Reid '93 


Nijm II 


Phenomenological 


67.8 


69.7 


71.1 


OPE (n = 5) 


69.8 


72.1 


74.0 
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A,(p) 




FIG. 11: (Color online) Coefficients of the expansion in Eq. (ID2[) . The red curve shows I = and 
the blue shows 1 = 2. 



using the OPE deuteron and the Avl8 initial state. 



1. First correction term: AO\ 



Consider the rightmost 5 in Eq. (j43 



AU i 7Z T77B 73 t -KNN\\V - k J 



-(p-k) 2 -/i 2 ^vivvvr y ;_ m ^ 



P 2 mi I 7 



x — 2mN - : / s-p, (Di) 

where /i 2 = 3m 2 r /4 and F^-atat is the form factor described in Appendix 3- The easiest way to 
evaluate the matrix element of this operator is to let the OPE act to the left on the deuteron 
in position space. The resulting expression is then transformed to momentum space. We can 
expand the fraction in the integrand of Eq. (1D1I) into spherical harmonics (taking pj = z), 



P 2 rriTr I 7 

2m N 2 ~ 2 



£4(p)Y,,o(p), (D2) 



and note that only the I = 0, 2 terms will contribute to s-wave production. The expansion 
coefficients are shown in Fig. (TTJ Clearly the I = 2 term is small and we neglect it here to 
avoid the extra algebra involved with a J = 2 operator (resulting in the A 2 reduced matrix 
elements in the notation of Ref. {sj). We find 

^1 = -34%. (D3) 
M 
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0.5 1.0 1.5 2.0 
FIG. 12: (Color online) Effect of the square root on the OPE a) tensor and b) central radial 
functions. 

2. Second O^g correction term: A0 2 



Next consider the term, 
A0 2 



o-i • (p - k)cr 2 ■ (k - p) f q (( , )2> 

-(p-k) 2 - " 2 



//- 



-fSW((p - k) 5 



v /(p - k) 2 + m 2 -m» 



x 



f^-^ + ^S-p. 
V2m 7V 2 27 H 



This term has a modified OPE, 



d 3 q iq . r <T 1 -qo- 2 -(-q) 2 , 
re i'nNNK^l J 



(27T) 



q 2 + /i 2 



y/q 2 + m 2 



(Ti ■ Vcr 2 • V 



47T 



12tt 



[^(O + o-i-tr^r)]. 



(D4) 

(D5) 

(D6) 
(D7) 



In Fig. [I2~lwe compare the functions T^(r) and Y^(r) to traditional OPE which has fi in place 
of the square root in Eq. (1D5j) . We use the Schrodinger equation to replace the p 2 /2m^ 
with V(r) and evaluate the matrix element in position space to find 

AM 2 



M 



+50%. 



(D8) 



3. Third 0% B correction term: AC 3 



Calculating the effects of the 5 in the denominator of the OPE is difficult to do exactly 
due to the combination of momenta that appear, 

1 



An <Ti-(p-k)<r 2 -(k-p) 2 2 . 

- (V(p-k) 2 + m 2 +5(p)J +m 2 /4 



(D9) 
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0.5 1.0 1.5 2.0 
FIG. 13: (Color online) Effect of the 5 on the OPE (a) tensor and (b) central radial functions. 

(recall that <5(p) = p 2 /2m^ — ^i-k/2 + 7/2). Instead we will evaluate it for fixed values of 5 
which represent the deviation of p away from p^, 



o + = o(pi + m^) = = 0.45 m„ 



S_ = S(pi- m„ 



2m N 



2m 



-0.32 m^. 



N 



M 

AA4 3 (*-; 
M 



+16% 
-32%. 



(D10) 
(Dll) 



The modified tensor and central functions and Y% are shown in Fig [T3j We define the 
correction as 

AM 3 (5) = M{5) -M(0). (D12) 

We find, 



(D13) 
(D14) 



4. Summary of corrections 

For the purposes of estimating the net result we take the average of the estimates in Sec. 
[D3] and find 



AM 



tot 



M 



-34% + 50% - 8% = +8%. 



(D15) 



We have successfully shown that the corrections to the first two-body DW amplitude are 
small, and actually increase the amplitude which is already twice as large as the traditional 
impulse approximation. 
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Appendix E: Effect of the 5 terms: O^-q diagram 



The second two-body DW amplitude's corrections are evaluated exactly as in the previous 
sub-sections and we just display the results here, 

AM i 



M 

AM 2 
M 

AM 3 (S+) AM 3 (5 



M M 



-35% (El) 
-39% (E2) 
= +3%, (E3) 



with the net result 



AMot 



( , - -35% - 39% + 3% = -71%. (E4) 
M 

We see that the corrections to the approximation in Eq. (14"6T) are fairly large, but this has a 
negligible effect because the amplitude is already small compared to the first two-body DW 
amplitude. 



Appendix F: Heavy meson exchange 



Consider the loop on the left-hand side of Fig. [TH which is obtained by using OPE for 
the left K in Eq. ( 133]) and a exchange (the dominant intermediate-range mechanism) for 
the right K . Note that this loop only differs from Fig. [6] in two ways: the meson-nucleon 
vertex (here we consider only scalar-isoscalar) and the meson mass. We use a typical set of 
parameters [33], g\j^ = 7.1 and m a = 550 MeV. 

The result of integrating over energy will proceed exactly as it did with the pion resulting 
in the two diagrams shown on the right-hand side of Fig. [HI To interpret the J^b* term 
(again, neglecting stretched box diagrams), we absorb the sigma exchange into the initial 
state and no new term is added. However, in the J^b term, after absorbing the pion exchange 
into the final state, we are left with a new operator. The amplitude for this operator can be 



obtained from that of Fig. 9(b) with the following change: 

3 



9a 



2 3 



~^ fir ~^ (fir) 2 



+ 1 



fir 



fi a r 



(Fl) 



where fi 2 = 3m 2 / 4 and fi 2 



(m w /2) 2 . We find, 



^2B 



-7.24, 



(F2) 



28 




FIG. 14: Impulse approximation with distorted waves: initial state heavy meson exchange. Solid 
lines represent nucleons, dashed lines pions, and the double solid line a a meson. Crosses represent 
propagators which are absent due to the partial fractions decomposition. 



which is larger in magnitude than the pionic A™> b (with the same sign) but smaller than 
A^' a (with the opposite sign). Since m a is relatively large, we can safely ignore the two 5 
corrections that are competing with Ui and only need to evaluate 

= -27%. (F3) 

M 

One natural question is whether the static a exchange already present in the initial-state 
wave function is a sufficient approximation for the contribution considered in this section. 
To answer this question, we can evaluate the traditional impulse approximation with 

mf = \pi,P2) + GV ff \*) h (F4) 

where here we employ a static a exchange that is present (at least effectively) in the wave 
function. Using this initial-state wave function, we calculate 

M° 1B = f (4>\0^)1, (F5) 

and find the reduced matrix element, 

A° 1B = -3.3. (F6) 

Thus we see that the a exchange in the traditional impulse approximation is an underesti- 
mate (in magnitude) of the true non-static exchange dictated by the loop integral. 

Of course there is no a in traditional B^PT, so this section is simply telling us that to 
achieve high accuracy it is indeed important to use more than just simple pion exchange 
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when forming the original box diagram. Such a calculation is beyond the scope of this work. 
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